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We present the results of a numerical investigation of charged-particle transport across a synthe- 
sized magnetic configuration composed of a constant homogeneous background field and a multiscale 
perturbation component simulating an effect of turbulence on the microscopic particle dynamics. 
Our main goal is to analyze the dispersion of ideal test particles faced to diverse conditions in the 
turbulent domain. Depending on the amplitude of the background field and the input test particle 
velocity, we observe distinct transport regimes ranging from subdiffusion of guiding centers in the 
limit of Hamiltonian dynamics to random walks on a percolating fractal array and further to nearly 
diffusive behavior of the mean-square particle displacement versus time. In all cases, we find com- 
plex microscopic structure of the particle motion revealing long-time rests and trapping phenomena, 
sporadically interrupted by the phases of active cross-field propagation reminiscent of Levy-walk 
statistics. These complex features persist even when the particle dispersion is diffusive. An in- 
terpretation of the results obtained is proposed in connection with the fractional kinetics paradigm 
extending the microscopic properties of transport far beyond the conventional picture of a Brownian 
random motion. A calculation of the transport exponent for random walks on a fractal lattice is 
advocated from topological arguments. An intriguing indication of the topological approach is a 
gap in the transport exponent separating Hamiltonian-like and fractal random walk-like dynamics, 
supported through the simulation. 
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I. INTRODUCTION 

Anomalous transport phenomena in complex nonlinear 
dynamical systems can often be associated with the effect 
of turbulence on the microscopic particle kinetics. The 
multiscale interaction of nonlinear relaxation processes 
and self-organization mechanisms operating in turbulent 
media customarily results in a transition to a nonequilib- 
rium (quasi)stationary state [1] dominated by long-range 
correlations in space and time [2]. In statistically ho- 
mogeneous, isotropic systems, turbulent correlations at 
the (quasi)stationary state support self-similar transport 
regimes that are far beyond the conventional Gaussian 
diffusion. As a rule, such regimes explicitly involve the 
impacts of memory, intermittency, and nonlocality; ex- 
amples are Levy-type processes [3-6] and fractal time 
random walks (FTRW's) [3,7,8], dating back to Mandel- 
brot and Van Ness fractional Brownian motion [9,10]. 

Levy processes incorporate bursty dynamics with al- 
most instant multiscale jumps often referred to as 
"flights" [4]. The defining feature of Levy flights is a 
heavy tailed power-law jump length distribution leading 
to infinite second moments of the corresponding prob- 
ability density. The problem of infinite moments may 
be avoided by replacing Levy flights with Levy walks 
through a spatiotemporal coupling posing inertia restric- 
tions over the bursty events [11,12]. Levy walks pro- 
vide suitable probabilistic model of superdiffusive trans- 
port [13] corresponding to faster-than-linear growth of 
the mean-square particle displacement versus time. In 
contrast to Levy flights, FTRW's are continuous random 
walk processes without identifiable jumps. Fractal time 
random walk statistics relies on a power-law distribution 
of waiting times between consequent steps of the mo- 
tion [14]. FTRW's account for long-time particle rests in 
the turbulent medium and customarily result in slower- 
than-lincar evolution of the tracer dispersion. A compre- 
hension of the essential role played by FTRW and Levy 
statistics in the microscopic description of turbulence has 
led to a formulation of the fractional — "strange" — ki- 
netics [15] reviewed in Refs. [16,17]. 

Near a (quasi)stationary state, the competition be- 
tween fractal time and Levy walk regimes produces a 
self-similar trajectory with the variance which grows with 
time raised to a power between and 2: 

(r 2 (t)) - 2Vt», 0<^<2. (1) 

The quantity V in Eq. (1) has the sense of a general- 
ized transport coefficient and is measured in units cm 2 
s~ M . One distinguishes between subdiffusive [p, < 1) and 
superdiffusive (ji > 1) processes: The former are domi- 
nated by long-time rests, the latter, by "active" migra- 
tion regimes such as Levy walks. The property of self- 
similarity appears in fractal structure of the trajectories 
over a broad range of spatial scales. The key parameter 
describing the topology of the particle track in the real 
space is the fractal dimension of the motion, d w > 1 [18]. 



The value of d w is related to the transport exponent /i in 
Eq. (1) via 

fi = 2/d w , d w > 1. (2) 

Equation (2) holds both FTRW's and flights [16]. The 
transport regimes included in Eqs. (1) and (2) range 
from ballistic motion (d w = l,p = 2) to the confine- 
ment (d w = oo,fi = 0). Diffusive processes correspond 
to d w — 2 and fi = 1. A realization of a trajectory hav- 
ing d w — 2 is provided by the conventional Brownian 
random walk model [10] showing uncorrelated behavior 
on all time and spatial scales [16]. Note that the Brown- 
ian random motion on a plane covers densely everywhere 
the ambient two-dimensional space as time t — > oo [19]. 
In case of subdiffusion (p < 1), the trajectory covers the 
plane with excess density when compared to the random 
Brownian counterpart, meaning d w > 2 in Eq. (2). The 
phenomenon is generally agreed to be due to multiscale 
memory effects which make the tracer return more often 
to the points already visited to time t. In contrast to 
subdiffusion, superdiffusive regimes (/i > 1) suppress the 
returns, enabling "fast" tracks with the "low" dimension 
d w < 2. An inherent drive for superdiffusion is often as- 
sociated with spatiotemporal nonlocality [16], though its 
detailed comprehension is far from being complete [14]. 

In many cases, the subdiffusive (/i < 1) behavior in Eq. 
(1) is supported by a self-organized concentration of tur- 
bulent transport processes on low-dimensional (fractal) 
arrays such as percolating lattices. [The term "percolat- 
ing" is synonymous with "infinite connected," meaning a 
structure which stretches to arbitrary long (macroscopic) 
scales. In what follows, "connected" is understood in a 
somewhat restrictive sense "path-connected:" An arbi- 
trary point found on a path-connected fractal set can 
propagate to another one along a continuous trajectory 
which lies everywhere on this set. The property of path- 
connectedness makes it possible to consider the fractal 
distribution as a single topological object.] Given a self- 
similar percolating fractal lattice with the Hausdorff di- 
mension df > 1, one finds the mean-square displacement 
of the walker after time t to be [18] 

(r 2 (t)) ~2Vt 2 /<- 2+e \ (3) 

The dispersion in Eq. (3) corresponds to the exponent 
/i = 2/(2 + 6) and the fractal dimension of the motion 
d w = 2 + 9. The quantity 9 is the index of connectivity 
of the fractal, which simultaneously appears in the prob- 
ability to return to the starting point p{t) oc £~ d //( 2 + 6 ') 
[20]. The index of connectivity describes intrinsic topo- 
logical features of fractal objects [21] and observes re- 
markable invariance properties [22]. For path-connected 
fractal distributions, the value of 9 > (as opposed to 
disconnected fractals having 9 < 0) [22]. In the limit of 
Euclidean (nonfractal) geometry, the index of connectiv- 
ity 9^0, leading to the conventional, Einstein relation 
(r 2 (t)) ~ 2Vt for the particle diffusion. Together with 
the Hausdorff dimension df, the index of connectivity 
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9 defines the so-called spectral (or fracton) dimension 
d s = 2d//(2 + 6) = 2d f /d w = fj,d f [23]. The value of d s 
determines the effective (fractional) number of degrees 
of freedom on a fractal geometry [24]. This fractional 
number is manifest in the scaling p(t) oc t~ ds / 2 , associ- 
ated with the return probability p(t) in a fractional "Eu- 
clidean" space enabling d s "orthogonal" directions [22]. 
In view of > 0, the value of the spectral dimension 
cannot exceed its Hausdorff counterpart as soon as the 
property of path-connectedness applies, d s < df. 

Physical realizations of strange transport processes 
(both subdiffusive and superdiffusivc) encompass fluid 
[25] and electrostatic drift-wave turbulence [26,27], 
chaotic flows generated by multipoint vortices [28], self- 
avoiding random walks [29], disordered solid materials 
[30], astrophysical [31] and laboratory [32] plasma, in- 
cluding highly intermittent, self-similar regimes mea- 
sured in the edge and scrape-off layer region of fusion 
devices [33]. In cosmic electrodynamics, unconventional 
statistical properties of transport bridging to the strange 
kinetics paradigm have been speculated for the solar pho- 
tosphere [34], the Earth's dayside magnetopause [35,36], 
and the Earth's distant magnetotail [37,38]. 

In this paper, we analyze self-similar transport pro- 
cesses in a synthesized magnetic configuration composed 
of a constant homogeneous background field and a multi- 
scale perturbation component simulating a marginal im- 
pact of turbulence on the behavior of test-particles. The 
physics included in the consideration below ranges from 
laboratory (fusion) realizations to typical conditions in 
space plasmas, specific to charged particle (ion) dynam- 
ics in turbulent magnetotail-like current sheets. Our 
study includes both numerical and theoretical counter- 
parts. The details of the numerical magnetic field model 
along with the microscopic equations of motion are given 
in Sec. II. The basic transport regimes a priori specu- 
lated for the turbulent system concerned are addressed 
in Sec. III. Numerical simulation results for the particles 
faced to diverse conditions in the turbulent domain are 
presented in Sec. IV. We end with a summary in Sec. V. 

II. MAGNETIC FIELD MODEL 

Setting (x, y, z) to be Cartesian coordinates in a three- 
dimensional Euclidean space, we define the magnetic field 
configuration through B x = 0, B y = 0, and B z = B, 
where 

B = B + 6B(r), v=(x,y,0), (4) 

B is a constant homogeneous background field, SB(r) is 
a magnetostatic perturbation generated across the hori- 
zontal (xy) plane, and B = (0, 0, B) is the magnetic field 
vector looking everywhere in the direction normal to the 
plane, z. The existence of a zero- frequency magneto- 
static perturbation mode in a magnetized plasma was 



shown in Ref. [39] . This mode is reminiscent of the zero- 
frequency electrostatic convective-cell (vortex) mode [40] 
in two-dimensional plasma. In particular, the two modes 
share scaling properties of the cross-field transport coeffi- 
cient above the classical collisional limit [39] . The behav- 
ior of extremely low-frequency perturbation modes in an 
unmagnetized plasma have been discussed by Ginzburg 
and Ruhadze [41]. The perturbation component in Eq. 
(4) is assumed in the form 

SB(r) = <5£>k exp[ik • r + itfk], (5) 

k 

where SB^ is the Fourier amplitude of the magnetostatic 
mode with wave vector k = (k x , k y , 0), and {<^k} are ran- 
dom phases posing — in the limit of statistical descrip- 
tion — the reflection symmetry 8B(r) <-> —6B(r). This 
mirrors the inherent parity properties of the multiscalc 
magnetic field turbulence in situ measured in the Earth's 
magnetotail current sheet [42]. Owing to the reflection 
symmetry, the separatrix B = Bq corresponding to the 
zero-set of the perturbation component SB(r) = con- 
tains a path-connected percolating subset [21] enabling 
field-line and magnetized charged-particle transport at 
the macroscopic scales [43]. The amplitudes SB^ intro- 
duced in the Fourier expansion in Eq. (4) are further 
defined through an algebraic function 

k ~ (fc 2 L 2 + l)(«+i)/4' W 

where k = |k| is the absolute value of the wave vector 
k, L is the side of the square, double periodic simula- 
tion box, A is the normalization constant, and a is the 
slope of the Fourier power-law energy density spectrum, 
P(k) <~ k~ a . Equations (5) and (6) define the so-called 
fractional Brownian surfaces [10], synthesized fractal ob- 
jects incorporating generic features of percolation in ir- 
regular media [43]. [Here we ignore nonlinear effects like 
mode coupling and build-up of correlations which may 
destroy the wave-like features in Eq. (5). The role of 
nonlinearity has recently been emphasized in Refs. [27] 
and [44] in connection with the formation of coherent vor- 
tical structures in low-/3 plasmas, dominating the trans- 
port on individual flux-surfaces. In this study, we focus 
on the impact of fractal behavior on the microscopic par- 
ticle dynamics, already present in the generic realization 
in Eqs. (5) and (6). ] 

In the numerical model, the components k x and k y are 
settled on the grid k x = 2im x /L and k y = 2im y /L, with 
integer valued n x and n y ; these satisfy n 2 ^ < n x +n y < 
n 2 lax . The k space is thereby a circular corona stretch- 
ing from fc min = 2irn min /L to fc max = 2irn maK /L. For the 
runs presented below, n m i n = 4 and n max = 80. Accord- 
ingly, the turbulence wavelengths A = 2ir/k range from 
Amin = L/80 to A max = L/4. The shortest wavelength 
Amin determines the finest scale of the inhomogeneities 
present, a <~ A m ; n . The longest wavelength A max , in its 
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turn, simulates the macroscopic turbulence correlation 
length, £ ~ A max . The ratio £/a ~ 20. This mimics the 
properties of the turbulence in the Earth's stretched and 
thinned magnetotail near a marginal (quasi)stationary 
state [38]. The number of independent Fourier modes 
in Eq. (5) is, by order of magnitude, N <~ 7r(n^ ax — 
n 2 nin )/2 <~ 10 4 . A zero cut SB(r) = of the fractional 
Brownian surface in Eq. (5) is illustrated in Fig. 1. 

The equation of motion for a particle of mass m and 
charge q, migrating across the magnetic field B, reads 



m 



dv 



v x B 



(7) 



where v = dr/dt = (dx/dt,dy/dt,Q) is the velocity at 
point r = (x,y,0) at time t. Denote p <~ vmc/qB to be 
the particle Larmor radius in the magnetic field B. (Here 
v = |v| is the absolute value of the velocity vector v.) 
Let us further introduce dimensionless parameters I ~ 
p/a and b ~ 5B/B - The effective cross- field transport 
regime may be sensitive to t and b, as we now proceed to 
show. 



III. BASIC TRANSPORT REGIMES 

A. Hamiltonian limit 

Assume first I — > +0 and b <C 1, meaning magnetized 
particles everywhere in the (xy) plane. As is well known 
[45], the dynamics of magnetized particles is adiabatic, 
enabling one to rely on a guiding center — drift — approx- 
imation of the long-time (t 3> mc/qB ) particle motion. 
The guiding center drift velocity is given by [45] 



BxVB 



u 



2q Bo 



leading to a system of Hamiltonian equations 



dx 
~dt 



d<f>(x,y) 
dy 



dy 

dt 



d<S>(x,y) 
dx 



(8) 



(9) 



Here, &(x,y) = — (v 2 /2Bofl)SB(x,y) is time indepen- 
dent — stationary — Hamiltonian, and £1 <~ qB^/mc is 
a characteristic particle gyrofrequency. The Hamiltonian 
<S>(x,y) is considered as the zero frequency limit (u) — > 0) 
of a time varying, single frequency (oc e lu)t ) Hamilto- 
nian &(x,y, t). Starting from the periodic Hamiltonian 
$(x, y, t) and turning the perturbation frequency to zero, 
w — > 0, one may demonstrate [19] that the main contribu- 
tion to the cross-field transport comes from a thin layer 
surrounding the separatrix B = B 0} where the particle 
excursion periods along percolating isoenergetic contours 
resonate with the basic period of the field, 2tt/uj — > oo. 

In many ways, the Hamiltonian dynamics in Eqs. (8) 
and (9) resembles the convection of a magnetized plasma 



by Ex B drifts caused by zero- frequency electrostatic vor- 
tex modes across the externally applied confining mag- 
netic field [40] . The ubiquitous zero frequency — static — 
limit v — > extremizes the case of "strong" turbulence 
[46] , conventionally described by the so-called Kubo num- 
ber Q ~ u/toa 3> 1, where u — |u| is the absolute value 
of the vector u, and a/u is a typical migration time in the 
field. For Q ^> 1, the large-scale behavior of the turbu- 
lent transport coefficient in Eq. (1) can be summarized 
by the power-law dependence [47] 



(10) 



where the fractal dimension d w > 1 stands for the topol- 
ogy of the transport process in the real space. The factor 
w 2/d„ Den irid a 2 accommodates the strange (cx t 2 l d ™) 
evolution of the variance (r 2 (i)) versus time t. In case 
of a diffusion (d w = 2), the coefficient in Eq. (10) scales 
with Q as Q 2 / 3 . This contrasts the widely known Bohm 
scaling V oc Q 1 , which is traditionally affiliated with 
charged-particle diffusion in low-frequency wave fields, 
since the pioneering works in Refs. [39] and [40]. The 
deviation from the Bohm scaling mirrors the role of per- 
colation properties in the guiding center picture of the 
cross- field particle transport [48]. 

As the perturbation frequency vanishes, ui — > 0, the 
Kubo number in Eq. (10) diverges as an inverse, Q — > 
oo. The transport coefficient V cx u?l d ™ Qdw/(2d w -i) thcn 
goes to zero if the fractal dimension of the motion satisfies 
2/d w > d w /(2d w - 1), yielding 1 < d w < 2 + \/2. For 
d w — > 2 + \/2 w 3.41 [to be associated with the stationary 
Hamiltonian $(x,y) in Eq. (9)], the value of V saturates 
at 



V 



,V2„,2-\/2 



/2 



(11) 



and does not depend on oj. The corresponding transport 
law (r 2 (i)) oc f is a subdiffusion with the exponent 



= 2-V2w0.58 (t-»oo). 



(12) 



The subdiffusivc regime in Eq. (12) is dominated by 
long-time particle rests near the points of equilibrium, 
defined by dQ/dx — and d^/dy = 0. This behavior is 
asymptotic: t — > oo. In view of Eq. (12), the coefficient 
in Eq. (11) observes anomalous scaling with the magnetic 
perturbation b <C 1, i.e., V oc b 2 ^^ 2 . 

It is instructive to emphasize that the parameter p — 
2 — \[2 is the exact lower bound on the transport expo- 
nent for stochastic Hamiltonian systems with 1 ^ degrees 
of freedom [47] . This lower bound is determined by the 
extremely long particle rests near the points of equilib- 
rium, where the guiding center drift velocity vanishes. 
The prevalence of such long rests in the zero frequency 
limit suppresses the effect of topology of the separatrix 
on the cross-field transport rate. This appears in the fact 
that the value p = 2 — \[2 does not depent on the details 
of the geometry of the array B — B , nor on the way the 
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separatrix is folded in the ambient space. Such details 
come into play as soon as the magnetized (adiabatic) be- 
havior is relaxed. 



B. Random walks on a fractal array 

Near the points of equilibrium, the Hamiltonian ap- 
proximation in Eq. (9) is invalidated already for very 
small but finite values of the particle Larmor radius, 
£ ~ p/ a . Nonadiabatic effects are more pronounced in 
magnetic configurations with larger perturbation com- 
ponent, b <~ SB/B . The limitation b <C 1 assumed in 
Eqs. (8) and (9) may now be loosen to b < O, where O 
is a constant of the order of 1. Even if slight (£ <C 1), 
nonadiabaticity is important as it revives the particles at- 
tempting to rest at equilibrium. The phenomenon helps 
the charge carriers to keep up their mobility near the 
points of occasionally small magnetic gradients VB — ► 0. 
The cross-field particle migration is naturally enhanced 
in this case, implying \i > 0.58 in Eq. (1). 

The inclusion of slight nonadiabaticity allows the par- 
ticles to pass more readily through the points of equilib- 
rium and thereby wander along the separatrix B = Bq 
in an almost casual way. A suitable approach to the 
problem could be found within a class of random walk 
models [18] associated with the percolative geomety of 
the set B = S [21]. 

In the case of an isotropic, random-phased perturba- 
tion such as the fractional Brownian surface in Eqs. (5) 
and (6), the separatrix B = Bo forms a percolating array 
which observes, in addition to the path-connectedness, 
fractal properties in the range of scales between a ~ A m ; n 
and £ ~ A max (see Refs. [10] and [49]). Here, £ has the 
sense of the (upper) fractal correlation length which is 
known to diverge at criticality, £ — > oo [18]. For "phys- 
ical" fractals, the value of £ is of course finite (£ < oo), 
though far longer than a. The finitcness of £ poses an 
upper bound where the fractal geometry of percolation 
crosses over to a statistically homogeneous distribution, 
d s — > df — > 2. 

In two ambient dimensions, the percolation along the 
separatrix is critical, meaning a threshold at the level 
B = Bq. Crossing the threshold one changes the domain 
which contains the paths to infinity [43]. The formation 
of critical percolating fractal structures has recently been 
conjectured as a generic property of the (quasi)stationary 
states in complex nonlinear dynamical systems far from 
thermal equilibrium [22]. 

The criticality character supports universal behavior 
of the percolation transition, manifest in a variety of in- 
triguing properties such as independence of the type of 
the percolation problem and of the microscopic details of 
the lattice [43]. In what follows, we are interested in the 
universality of the spectral dimension d s [23], first conjec- 
tured by Alexander and Orbach [50] and later addressed 
in an improved form by Milovanov [51] who introduced 



the notion of the percolation constant C, a topological 
parameter incorporating the features of connectedness of 
fractal distributions. The percolation constant is defined 
as the smaller (between the two possible) root to the 
identity [51] 



T C/2 



r(c/2 + i) 



71", 



(13) 



where the symbol T denotes the Euler gamma function. 
The value of C, a transcendental parameter approxi- 
mately equal to 1.327 . . ., determines the least fractional 
number of degrees of freedom, enabling a particle to reach 
the point at infinity through a random walk process on 
a self-similar fractal geometry. A topological approach 
to the phenomenon of percolation leading to the identity 
in Eq. (13) is discussed in some detail in Ref. [22]. In 
terms of the percolation constant, the universality of the 
spectral dimension at the threshold is quantified by [51] 



d s = 2d//(2 + 0) =C w 1.327. 



(14) 



Once the number of degrees of freedom is known from Eq. 
(14), an evaluation of the exponent fi = 2/(2 + 0) = C/df 
in Eq. (3) is reduced to a calculation of the Hausdorff 
dimension df of the lattice on which the transport process 
concentrates. For critical percolation on two-dimensional 
fractal arrays, the value of df lies within [52] 



C < d f < S = In 8/ In 3 = 1.89. .. < 2, 



(15) 



where S = In 8/ In 3 is the Hausdorff dimension of the 
square Sierpinski carpet [53] , a celebrated topological ob- 
ject providing the universal embedding for paths on a 
plane [54]. Inequality (15) derives as a condition which 
brings together the property of path-connectedness and 
the threshold character, manifest in the fact that both C 
and S pose the restrictions on the Hausdorff dimension 
df. In view of \i = C/df we have, at the critical range, 



<M<1, 



1+n 



= C/S « 0.70. 



(16) 



(17) 



Note a gap A/z w 0.12 between /i m ; n = C/S w 0.70 and 
the Hamiltonian bound /j, = 2 — \/2 rts 0.58 in Eq. (12). 
Owing to the gap, a transition from Hamiltonian dynam- 
ics to random walks on a fractal pattern may have an 
abrupt character, when a slight increase in the particle 
Larmor radius I results in an almost sudden growth (from 
w 0.58 up to at least 0.70) of the transport exponent [i. 
This makes the Hamiltonian regime in Eq. (12) appre- 
ciably sensible to the particle Larmor radius t assumed 
in the simulation. 

The exact value of the Hausdorff dimension df, specific 
to the topology of fractional Brownian surfaces, has been 
calculated in Ref. [21] from the shape of fractal isoener- 
gctic contours in vicinity of the separatrix B = B : 
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df = 2C - 1 w 1.65. (18) 

The Hausdorff dimension in Eq. (18) falls well inside 
the range defined by inequality (15). Relation (18) cor- 
responds to the exponent of the Fourier energy density 
spectrum obeying 1 < a < 7 — 4C« 1.69 [21]. The latter 
guarantees enough energy at the small scales, validating 
self-similar behavior of the microscopic kinetic process. 
Combining Eqs. (3), (14) and (18), one arrives at 

H = C/(2C - 1) w 0.80. (19) 

The ensuing cross-field particle transport is a subdif- 
fusion, ^r 2 (i)^ oc t 080 . This is enhanced when com- 
pared to the subdiffusive behavior (r 2 (i)} oc t 58 in the 
limit of Hamiltonian dynamics quantified by Eqs. (8) 
and (9). The fractal dimension of the motion equals 
d w = 2/fj, = 2(2C -1)/C» 2.49 and is in fact lower than 
in the adiabatic case (d w = 2 + \[2 w 3.41). A derivation 
of the turbulent transport coefficient V for the particles 
having small but finite Larmor radius I <C 1 involves 
the issue of the branching dimension [19] and is intricate 
somewhat [35]. Here we note that V may not reveal a 
single scaling behavior, contrary to the Hamiltonian co- 
efficient in Eqs. (10) and (11). 

A subdiffusion consistent with the exponent in Eq. 
(19) has earlier been proposed for ions migrating across 
the turbulent magnetotail current sheet [38]. The esti- 
mate in Eq. (19) applies to the "anomalous" time scales 
a 2+e <t< £ 2+e for which the particles walk on the frac- 
tal. For longer times t > £ 2+e , the particles cover dis- 
tances generally exceeding the fractal correlation length, 
£. This implies a transition to a random walk process on 
an Euclidean (nonfractal) geometry. Hence 

fj, -» 1 (20) 

for t > £ 2+e . Consequently, the asymptotic (t — ► oo) be- 
havior of the cross-field transport process in the slightly 
nonadiabatic case (£ <C 1 , b < O) should be diffusive. 
The corresponding particle diffusion coefficient is of the 
order of 

V ~ au/2. (21) 

The occurrence of the asymptotic (t — > oo) diffusive 
regime already in the slightly nonadiabatic realization is 
in contrast to the magnetized particle dynamics relying 
on the Hamiltonian Eqs. (9). The latter is subdiffusive 
(with the exponent fi w 0.58) on all time scales t 3> I/O, 
up to t — > oo. Note that the diffusion cross-over time 
scale for random walks on a percolating fractal array di- 
verges as £ 2+e for £ — ► oo. 

C. Strong nonadiabaticity 

In the limit of strong (£ > 1) nonadiabaticity, the cross- 
field transport processes occupy wide stochastic domains, 



whose measure is comparable with that of the ambient 
Euclidean plane itself. The microscopic particle dynam- 
ics thereby evolves into random walks in two dimensions 
(d s — > df — > 2), meaning diffusive behavior fi — > 1 on the 
time scales t » a/ v. The cross-field diffusion coefficient 
can be evaluated as [cf. Eq. (21)] 

V ~ avi/2. (22) 

Remark that the condition £ > 1 is always satisfied near 
the zero-set SB(r) = for vanishing background com- 
ponent, B — > 0. The perturbation parameter diverges 
in this case, b — > oo, signifying a decaying role of the 
percolation properties. 

IV. NUMERICAL RESULTS 

In this section, we present a selection of numerical runs 
which illustrate the basic transport regimes operating in 
the magnetic configuration in Eqs. (4) — (6). Magnetic 
field models with both nonvanishing and zero background 
component are investigated. All lengths are normalized 
to the side of the simulation box, L, and all times, to the 
inverse of the particle gyrofrequency, f2. The running 
value of refers to the local magnetic field B which may 
appreciably deviate from B . The particle velocity v is 
measured in units £IL, with typical values ranging — for 
different realizations — from 2 x 10~ 4 to 5 x 10~ 4 . For 
each run, 5000 particles are injected in a random manner 
throughout the simulation box, L x L. Initial velocity 
vectors v are randomly directed. The equation of mo- 
tion (7) is solved numerically by means of a fifth-order 
Runge-Kutta integrator with adaptive step; the maxi- 
mum time step is 10~ 2 . Depending on the strength of 
the background field Bq and the exact value of the in- 
jection velocity v, the total integration time varies from 
3 x 10 6 to 1 x 10 7 . The accuracy of the computation is 
checked by various methods, including conservation en- 
ergy verification at the end of the run, showing relative 
errors of less than 10~ 4 . The evaluation of the exponent 
/j is based on running fits with a time window At ~ 10 6 , 
which is made gradually move along t. 

A. Zero background component 

A nearly diffusive regime is readily recovered for the 
zero background component Bq = (i.e., for the infi- 
nite perturbation parameter b — oo). The details of the 
simulation are as follows. The particles are injected at 
random with the velocities v = 5 x 10~ 4 . The integra- 
tion is performed up to 3 x 10 6 , enabling a rich statistics 
of tracking data. The value of the transport exponent 
computed over the whole set of the trajectories is found 
to be /i i=a 0.96 ± 0.05. This value incorporates the parti- 
cles trapped on closed isoenergetic contours far from the 
separatrix, B = B . Simultaneously, the dispersion of 
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those particles initialized in close vicinity of the zero-set 
SB(r) = and propagating to longest distances across 
the field is not distinguished from a diffusion-type pro- 
cess. Sample trajectories of such particles bringing a ba- 
sic contribution in the cross-field transport are illustrated 
in Figs. 2 and 3. We draw attention to the very com- 
plex microscopic structure of the particle motion, reveal- 
ing a chaotic alternation of long-time rest and sporadic 
transient periods, manifest in the highly irregular, inter- 
mittent way the trajectory is marked by the tracer. In 
this connection, the diffusion appears to be an intricate 
balance between the time intervals when the particles 
stay in traps, and the "green-light" regimes — reminis- 
cent of Levy walk statistics — when the charge carriers 
actively migrate across the medium. This intriguing pic- 
ture challenges the conventional Brownian random mo- 
tion paradigm customarily associated with diffusion. We 
interpret this unusual behavior as a "strange diffusion" 
[22] deriving from a parity between rests and walks, the 
two competing counterparts of the dynamics underlying 
the transport at the microscopic scales. 

B. Finite background component 

We turn next to a situation with a finite background 
field B . To support the effect of the particle Larmor 
radius, we set b ~ 2 in the simulation below. This choice 
is a suitable compromise between the opposite extremes 
of b <C 1 and b = oo. 



1. Random walks on a fractal array 

As in the regime with zero background component, we 
set the injection velocity v to 5 x 10~ 4 . Accordingly, the 
integration covers time scales up to 3 x I0 6 and repeats 
the details of the run specified in subsection A. The par- 
ticle dispersion versus time is plotted in Fig. 4, solid 
line. The fit yields fi w 0.84 ± 0.05. This value is clearly 
subdiffusive and complies with the exponent in Eq. (19). 
Consequently, we consider the estimate \i w 0.84 ± 0.05 
as an evidence for a concentration of the cross-field trans- 
port on a percolating fractal array owing to the departure 
of Bq from the zero limit. 

To help judge the result obtained, we launch another 
run with all the same integration and test particle pa- 
rameters, but with a filter in the injection scheme. As 
the filter is applied, the particles found "too far" from 
the percolation level B = B are discarded. Thus, we 
only integrate along the trajectories that start in a band 
B ± AS, with AB a fraction of B . We could thereby 
increase the portion of particles staying close to the per- 
colation level and appreciably speed up the computation. 
The behavior of the particle dispersion for this run is 
shown in Fig. 4, dashed line. While larger distances are 
achieved, the exponent \i remains practically unchanged, 



H w 0.86 ± 0.05. This observation proves that the trans- 
port, as a matter of fact, is mostly due to the particles 
migrating in close vicinity of the percolation level. 

A realization of the percolation-associated particle tra- 
jectory in the presence of the background component is 
developed in Fig. 5. The tracer features a tendency to- 
ward adiabatic dynamics, manifest in the periods of a 
drift-like migration along the separatrix. Near the sad- 
dles of the zero-set SB(r) = (see Fig. 1), the drift-like 
regimes are interrupted by almost unmagnetized, mean- 
dering motion in a close-to-separatrix layer owing to the 
effect of finite Larmor radius. The observation of mean- 
dering is important as it supports the chance for casuality 
enabling to consider the cross-field transport in connec- 
tion with random walk processes on percolative fractal 
geometry [18]. At the microscopic scales, the meander- 
ing motion bears signatures reminiscent of the intermit- 
tent dynamics found for B = (see Figs. 2 and 3). 
The drift, in its turn, exhibits inhomogeneous behavior 
with alternating phases of fast and slow progress. These 
competing phases of fast and slow migration acquire the 
role of the walk-and-rest statistics as the adiabatic limit 
is approached. We emphasize that the "walks" are in- 
herently present in the motion, even though the ensuing 
transport appears to be subdiffusive. In this regard, the 
parity between walks and rests previously speculated for 
the strange diffusion is now shifted toward rests. 

2. Back to diffusion 

With increasing velocity v, the mean-square displace- 
ment of the tracer tends to a linear (diffusive) form start- 
ing from longer integration times. For shorter times, the 
dispersion may still be subdiffusive. For instance, turn- 
ing v from 5 x 10~ 4 to 1 x 10~ 3 , we locate the transport 
exponent within fi w 0.95 ± 0.05 for the long integration 
time 1 x 10 7 . As the integration time is reduced to 3 x 10 6 , 
the fractal random walk-like dispersion fj, w 0.84 ± 0.05 
is recovered. Runs with larger velocities reveal nearly 
diffusive behavior already for the short times 3 x 10 6 . 
These results demonstrate the effect of the finite fractal 
correlation length, £, posed by the basic periodicity of 
the simulation box. We emphasize that the periodic ex- 
tension of the simulation box assumed in the simulation 
truncates the percolative fractal geometry of the field be- 
yond £ ~ A max . 

3. Hamiltonian limit 

To approach the adiabatic regime in Eq. (9), we turn 
the injection velocity v to 2 x 10~ 4 . Simultaneously, we 
maintain the background component Bq at the level cor- 
responding to b ~ 2 in order to achieve a clearer compar- 
ison with the previous runs. A longer integration time 
1 x 10 7 matching the proposed decrease in v is settled. 
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The particle dispersion as a function of time is summa- 
rized in Fig. 6. 

As the velocity is set to 2 x 1CP 4 , a pronounced subd- 
iffusive behavior is observed through the integration pe- 
riod. The fit yields fi w 0.67 ±0.05, which is appreciably 
smaller than the values allowed for a random walk pro- 
cess on a fractal geometry [see Eqs. (17) and (19)]. The 
estimate of /i 0.67±0.05 is reminiscent of the exponent 
ji w 0.58 in Eq. (12), corresponding to the limit of Hamil- 
tonian dynamics. Some deviation between the numerical 
result /i w 0.67 ± 0.05 and the extreme of /j, w 0.58 could 
be due to the excessively large magnetic perturbation pa- 
rameter, b ^ 2. This deviation may further be shown to 
decrease with increasing background field, Bo. 

Turning, gradually, the injection velocity v from 5 x 
10~ 4 to 2 x 10~ 4 , we find an unstable transitional behav- 
ior when the exponent \x fluctuates between Hamiltonian- 
like and fractal random walk-like values. This poses an 
uncertainity in /i for some intermediate values of v. To 
the major extent, the uncertainities are bypassed as v 
gets sufficiently close to 2 x 10~ 4 . We associate the ob- 
served unstable regimes with the existence of the gap 
deriving from Eqs. (12) and (17). On the contrary, no 
specific unstable domain is found for v varying between 
5 x 10~ 4 and 1 x 10~ 3 , showing a continuous turnover to 
the diffusive transport. 

A sample trajectory corresponding to the "cold" par- 
ticles with velocities v = 2 x 10~ 4 is plotted in Fig. 7. 
The motion is drift-like everywhere along the isoenergetic 
contour. No unmagnetized, meandering effects are recog- 
nized, meaning a saturation of the transport process at 
the adiabatic regime. The drift, nevertheless, reveals in- 
homogeneities resembling the "walk-and-rest" statistics, 
manifest in the phases of "fast" and "slow" cross-field 
propagation. These signatures, already noticeable for 
larger injection velocity 5 x 10~ 4 , now entirely dominate 
the dynamics. 

V. SUMMARY AND CONCLUSIONS 

We have analyzed the dispersion of charged particles in 
a synthesized magnetic field configuration composed of a 
constant homogeneous background field and a multiscale 
perturbation component simulating the effect of turbu- 
lence on the test-particle dynamics. Already for such 
a simple magnetic field model, we observe complex mi- 
croscopic behavior of the charge carriers, demonstrating 
chaotic alternation of periods of rest and active cross- 
field migration regimes reminiscent of Levy-walk statis- 
tics. The integral process looks like a competition be- 
tween rests and walks governing the transport properties 
on the microscopic scales. 

In absence of the background component, the cross- 
field transport is almost diffusive, showing linear growth 
of the mean-square particle displacement versus time. At 
the microscopic level, this regime profits from an intri- 



cate compromise between the time intervals the parti- 
cles stay in traps and the periods they actively propa- 
gate through the medium. Such complex features are in 
contrast with the conventional Brownian random motion 
paradigm, often associated with diffusion. With increas- 
ing background field, the equilibrium between rests and 
walks — the two constituents of the motion — shifts to- 
ward longer rests, leading to a sublinear (subdiffusive) 
behavior of the particle dispersion. We emphasize that 
the walks are inherently present in the microscopic pic- 
ture of the dynamics, though the contribution they bring 
to the dispersion may be suppressed by the effect of rests. 
Two distinct subdiffusive regimes have further been rec- 
ognized, depending on the injection velocity v (i.e., on 
the particle Larmor radius, I ~ vmc/aqB). 

In the numerical model, a suitable range for the veloc- 
ity v is allocated — in dimensionless units — from 2x 10~ 4 
to 5 x 10~ 4 . For the larger values of v <~ 5 x 10" 4 en- 
abling moderate nonadiabaticity of the charge carriers, 
the cross-field transport is concentrated on percolating 
fractal arrays associated with the zero-set of the per- 
turbation component, 5B(r) = 0. In two ambient di- 
mensions, the percolation along the zero-set is critical, 
meaning a threshold at the level B — Bq. The value of 
the transport exponent as predicted by the percolation- 
based model equals /U w 0.80. This result — deriving 
from the fundamental properties of universality of the 
percolation transition — is reproduced in the simulation, 
[i w 0.84 ± 0.05. A characteristic feature of the percola- 
tion regime is the presence of chaotic meandering motion 
in a thin layer enveloping the separatrix B = Bq, mixed 
with gyromotion and trapping. 

As the velocity is turned to v = 2 x 10~ 4 , the trans- 
port exponent drops to ji rj 0.67 ± 0.05. This behavior 
is reminiscent of the static limit fi rj 0.58 of stochastic 
Hamiltonian systems with l| degrees of freedom. The 
tendency toward Hamiltonian dynamics reflects the pro- 
posed decrease in v in the presence of the background 
field. Some deviation between the numerical estimate of 
/i Rj 0.67 ± 0.05 and the extreme ofps; 0.58 may further 
be shown to be due to the excessive effect of the mag- 
netic perturbation introduced in the simulation. This 
effect decays with increasing background field. 

A main result of the analytical investigation — per- 
formed in parallel with the numerical study — is that 
the regime of Hamiltonian dynamics is separated from 
random walk processes on a percolative fractal geometry 
by a gap in the transport exponent, A^i rts 0.12. In the 
numerical simulation, the gap appears in an uncertain 
behavior of the particle dispersion for some intermediate 
values of 2 x 10~ 4 <d<5x 10~ 4 . Within this interval, 
the exponent \i fluctuates between Hamiltonian-like (« 
0.67 ±0.05) and fractal random walk-like (rj 0.84 ±0.05) 
dispersion. The unccrtanities are generally bypassed as 
v gets sufficiently close to 2 x 10~ 4 . 

As the Hamiltonian regime is approached, the mean- 
dering counterpart of the particle dynamics vanishes. In 
this limit, the cross-field transport is dominated by long- 
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time particle rests near the points of occasionally small 
magnetic gradients. Posing slight nonadiabaticity readily 
calls for the meandering to come into play. The trans- 
port exponent \i then promptly returns to the fractal ran- 
dom walk-like dispersion. With increasing nonadiabatic- 
ity the asymptotic diffusive regime is recovered starting 
from longer integration times. This result incorporates 
the effect of the finite fractal correlation length, posed by 
the basic periodicity of the simulation box. 
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FIG. 1. A zero cut SB(r) = of the fractional Brownian 
surface in Eqs. (4) and (5), corresponding to the exponent 
of the Fourier energy density spectrum a = 3/2. Owing to 
the sign parity SB(r) «-» —SB(r), the cut 8B(r) = contains 
a percolating isoenergetic contour enabling charged-particle 
transport across the field already in the limit of adiabatic 
dynamics. 

FIG. 2. A track of a test-particle injected in close vicin- 
ity of the zero-set SB(r) — and propagating to longest 
distances present. The background component is zero ev- 
erywhere, Bo = 0, corresponding to the infinite value of the 
perturbation parameter, b = oo. The injection velocity equals 
» = 5x 10~ 4 . Dimensionless units. 

FIG. 3. A typical test-particle trajectory underlying the 
"strange diffusion." All the same injection and magnetic re- 
alization parameters as in Fig. 2. Dimensionless units. 

FIG. 4. The particle dispersion as a function of time for the 
finite perturbation parameter, b ~ 2. The injection velocity 
set to v = 5 x 10~ 4 . The integration covers time scales up to 
3x 10 6 . Both runs with (dashed line) and without (solid line) 
the filter are shown. In absence of the filter, the particles 
are injected everywhere at random. The filter discards the 
particles found "too far" from the percolation level. 

FIG. 5. A track of a particle faced to random walk pro- 
cess on a percolative fractal geometry. Note the characteris- 
tic meandering motion mixed with inhomogeneous drift-like 
migration along the separatrix. 

FIG. 6. The particle dispersion versus time for two different 
values of the injection velocity, « = 5x 10~ 4 (solid line) and 
v — 2 x 10~ 4 (dashed line). Finite perturbation parameter, 
b ~ 2. The integration is performed up to 1 x 10 7 . 

FIG. 7. A track of a nearly magnetized particle migrating 
along the percolating isoenergetic contour. In contrast to the 
trajectory in Fig. 5, no meandering-like motion is present. 
Dimensionless units. 
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